In major dysmood disorder, physiosomatic, chronic fatigue and fibromyalgia symptoms are driven by immune activation and increased immune-associated neurotoxicity

Major depressive disorder (MDD) is accompanied by activated neuro-immune pathways, increased physiosomatic and chronic fatigue-fibromyalgia (FF) symptoms. The most severe MDD phenotype, namely major dysmood disorder (MDMD), is associated with adverse childhood experiences (ACEs) and negative life events (NLEs) which induce cytokines/chemokines/growth factors. To delineate the impact of ACE + NLEs on physiosomatic and FF symptoms in first episode (FE)-MDMD, and examine whether these effects are mediated by immune profiles. ACEs, NLEs, physiosomatic and FF symptoms, and 48 cytokines/chemokines/growth factors were measured in 64 FE-MDMD patients and 32 normal controls. Physiosomatic, FF and gastro-intestinal symptoms belong to the same factor as depression, anxiety, melancholia, and insomnia. The first factor extracted from these seven domains is labeled the physio-affective phenome of depression. A part (59.0%) of the variance in physiosomatic symptoms is explained by the independent effects of interleukin (IL)-16 and IL-8 (positively), CCL3 and IL-1 receptor antagonist (inversely correlated). A part (46.5%) of the variance in physiosomatic (59.0%) symptoms is explained by the independent effects of interleukin (IL)-16, TNF-related apoptosis-inducing ligand (TRAIL) (positively) and combined activities of negative immunoregulatory cytokines (inversely associated). Partial least squares analysis shows that ACE + NLEs exert a substantial influence on the physio-affective phenome which are partly mediated by an immune network composed of interleukin-16, CCL27, TRAIL, macrophage-colony stimulating factor, and stem cell growth factor. The physiosomatic and FF symptoms of FE-MDMD are partly caused by immune-associated neurotoxicity due to T helper (Th)-1 polarization and M1 macrophage activation and relative lowered compensatory immunoregulatory protection.

The specific mechanisms through which ACEs or ACE + NLEs could impact the physiosomatic and FF symptoms of MDD/MDMD remain uncertain.It is unclear whether these effects are mediated by the activation of the IRS, specific immune phenotypes like Th-1 or M1 phenotypes, or specific cytokines, chemokines, and growth factors such as IL-6, IL-1β, and TNF-α.Therefore, the precise hypotheses of this study are as follows: the impact of ACEs and NLEs on physiosomatic and FF symptoms is influenced by increased activity of the IRS and Immu-NT, specifically through the release of certain M1 cytokines such as IL-6, IL-1β, and TNF-α.

Participants
In accordance with the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) 40 , a senior psychiatrist at Alhakiem Hospital in Najaf, Iraq, identified and recruited 64 individuals with FE-MDD and 32 healthy controls.Moreover, all patients were diagnosed with MDD, single episode, indicating they experience their first MDD episode, in the current study labeled as FE-MDMD.All patients diagnosed with MDD were observed to be in the acute phase of the disease, with no indications of complete or partial remission.It should be noted that all FE-MDD patients met the diagnostic criteria for MDMD as detailed by Maes, Vasupanrajit, et al.  (2023)  16 .The diagnosis of MDMD versus SDMD is made using machine learning techniques (support vector machine) and as such there are no easy case definitions that can be applied.Nevertheless, MDMD is the more severe subtype in terms of higher ROI and/or severity of depression and anxiety.For example, when both the HAMD and the HAMA exceed twenty-two points, MDMD may be diagnosed even in patients with a first episode.
In addition, the senior psychiatrist recruited healthy controls from the same geographical area to serve as the control group.This group included hospital staff and patients' acquaintances.All participants were enlisted between October 2021 and March 2022.Exclusions were made for those with chronic liver or kidney disease, women who were pregnant or nursing, and those with a history of multiple sclerosis, Parkinson's disease, stroke, or Alzheimer's disease.Similarly, psoriasis, rheumatoid arthritis, inflammatory bowel disease, cancer, type 1 diabetes, and scleroderma exhibited the same pattern.In addition, we excluded participants who had experienced an acute COVID-19 infection, severe or critical COVID-19 disease, Long COVID, or a COVID-19 infection within the previous six months.Our study did not include individuals taking immunosuppressive or immunomodulatory medications, nor did it include those taking therapeutic doses of antioxidants or omega-3 supplements (last three months).
The process of meticulously selecting the control group was accorded significant consideration.Excluded from the control group were participants with a documented lifetime history of clinical depression or dysthymia, a family history of depression, mania, psychosis, or substance use disorder, or a history of suicide.Dysthymia (excluding cases of double depression), schizophrenia, schizoaffective disorder, bipolar disorder, autism spectrum disorders, substance use disorders (excluding nicotine dependence), post-traumatic stress disorder, psychoorganic disorders, generalized anxiety disorder, and obsessive-compulsive disorder were excluded from the study.
Before participating in the study, all participants, or, if applicable, their parents or legal custodians, provided informed consent in writing.Document No. 18/2021 indicates that the Ethics Committee of the College of Medical Technology at the Islamic University of Najaf in Iraq has approved this investigation.The research was conducted in accordance with both Iraqi and international ethical and privacy regulations.In prominent documents such as the Declaration of Helsinki by the World Medical Association, the Belmont Report, the CIOMS Guideline, and the International Conference on Harmonization of Good Clinical Practice, a variety of non-exhaustive principles are outlined.Our organization's institutional review board (IRB) is committed to upholding the highest standards of quality, ensuring precise compliance with the International Guideline for the Conduct of Safe Human Research (ICH-GCP).

Clinical assessments
The physical, mental, and behavioral health of the participants was evaluated by a senior psychiatrist who used a systematic interview and standardized procedures.The senior psychiatrist collected sociodemographic, clinical, and psychological information through semi-structured interviews.The same professional utilized the HAMD and the HAMA to evaluate the severity of depression and anxiety, respectively.In the current study, all physiosomatic symptoms of the HAMD and HAMA were excluded to compute pure depression and pure anxiety scores, respectively.The former concept was conceptualized as a collection of symptoms including depressed mood, feelings of guilt, suicidal ideation, and decreased interest.The score for pure anxiety was determined by adding the scores for anxious mood, tension, fears, and anxious behavior observed during the interview.The physiosomatic symptom score was calculated as a z unit-based composite score based on the sum of the z scores of HAMD and HAMA physiosomatic symptoms, namely anxiety somatic, somatic GIS, and genitourinary symptoms, hypochondriasis, somatic sensory, cardiovascular, GIS, genitourinary, autonomic, and respiratory symptoms.The pure FF symptoms were calculated as the sum of the FF items after exclusion of non-physiosomatic symptoms: muscle pain, muscle tension, fatigue, autonomic, gastro-intestinal symptoms, headache, and a flu-like malaise (all FF scale items after exclusion of all non-physiosomatic symptoms, including cognitive deficits, and sadness).GIS symptoms were conceptualized as a z unit-based composite score based on symptoms of the HAMA and FF scores, namely somatic, and GIS symptoms.Insomnia was conceptualized as a z value-based composite, namely sum of insomnia (HAMA item), insomnia early + insomnia middle + insomnia late (HAMD items) + sleep disorders (FF).Melancholia was conceptualized as a composite based on the z scores of insomnia late, psychomotor retardation, psychomotor agitation, loss of weight, and diurnal variation.This study utilized two items of the Columbia Suicide Severity Rating Scale (C-SSRS) to evaluate suicidal behaviors: the quantification of suicide attempts within the previous year and the assessment of the frequency of suicidal ideation within the previous three months 41 .A composite score of suicidal behaviors (SB) was calculated by adding the weighted z scores for the HAMD suicide item, the number of suicide attempts, and the frequency of suicidal ideation.
The ACEs Questionnaire 42 was used to measure the extent of adverse childhood experiences.There are a total of twenty-eight items on the scale, which comprise the scoring of ten distinct domains.These include (1) mental trauma, (2) physical trauma, (3) sexual abuse, (4) mental neglect, (5) physical neglect, (6) witnessing domestic violence involving the mother, (7) presence of a family member with drug abuse issues, (8) presence of a family member with depression or mental illness, (9) experiencing the loss of a parent due to separation, death, or divorce, and (10) having an incarcerated family member.Various ACE scores were calculated, including the total sum of all ACEs (termed total ACE), the sum of physical trauma, mental neglect, and family member with substance abuse (termed ACE247), and the sum of ACE247 and a family member in prison (termed ACE24710).In addition, the NLEs scale was used to assess the incidence of NLEs in the preceding year 43 .For the purposes of this investigation, the following factors were considered: a serious accident, the death of a family member or close friend, divorce or separation, unemployment, job loss, alcohol-related issues, drug-related issues, witnessing physical altercations or assaults, experiences of abuse or violent crime, encounters with law enforcement difficulties, problem gambling, familial incarceration, overcrowding in the household, and instances of discrimination.Consequently, we performed calculations to ascertain the combined effects of ACE and the occurrence of one or more NLEs, which we referred to as ACE + NLEs.Body mass index (BMI) was calculated by dividing the participants' weight in kilograms by the square of their height in meters.Using the diagnostic criteria specified in the fifth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5), tobacco use disorder was evaluated.

Biochemical assays
In this investigation, a disposable syringe and serum tubes were used to obtain a 5 mL venous blood sample from each participant while they were fasting.Blood samples were collected between 8:00 and 11:00 in the morning.The serum and blood cells were effectively separated following centrifugation at 35,000 revolutions per minute (rpm).The serum was then transferred into small Eppendorf containers and stored at − 80 °C until its use in the biomarker assays was deemed necessary.Bio-Plex Multiplex Immunoassay kits (Bio-Rad Laboratories Inc., Hercules, USA) were used to measure the levels of forty-eight cytokines/chemokines/growth factors in the serum of all participants.For measuring the concentrations of these proteins in blood serum, we employed a fluorescence-based method.The researchers measured the fluorescence intensity (FI) of each protein and used the values obtained after subtracting the blank value 24,44 .The Electronic Supplementary File (ESF) comprises Table 1, which provides an exhaustive listing of the analytes identified by our research.In addition, each analyte's respective gene ID and alternate names are listed in this table.The coefficients of variation (CV) for all analytes in the assay were less than 11.0%.To determine the concentrations of the analytes, we employed the manufacturersupplied standards.Subsequently, the out-of-range (OOR) concentration rate was calculated, i.e., the proportion of concentrations that fell below the minimum detectable level.Cytokine/chemokine/growth factor levels were excluded from the statistical analysis when more than 80% of the concentrations were below the lowest OOR value.Consequently, IFN-α2, IL-3, IL-7, and IL-12p40 were excluded from our study.Prevalences (measurable or not) were assigned to analytes whose levels exceeded the OOR in 20% to 40% of the cases.These variables were incorporated into the analysis as dummy variables (0 = not measurable, 1 = measurable).These prevalence rates were taken into account in the calculation of our composite immune profile scores.As a result, diverse immune profiles were computed, as described in Table 2 of the ESF.These profiles covered all analytes with the exception of those with detectable concentrations below 20% 24,44,45 .Table 2 of the ESF lists the variables used to generate M1, M2, M1/M2 (z M1-z M2), Th1, Th2, Th1/Th2 (z Th-1-z Th-2), Th-17, IRS, CIRS, IRS/CIRS (z IRS-z CIRS), and Immu-NT profiles.In addition, the individual cytokines, chemokines, and growth factors were analyzed to ascertain whether there were significant differences in immune profiles between the study groups.

Statistics
This study's statistical analyses were performed with IBM SPSS 29, Windows version.To compare continuous variables between study groups, statistical tests such as analysis of variance (ANOVA) and the Kruskal-Wallis test were utilized.Multiple post-hoc comparisons (Fisher's exact probability test) were performed to examine differences between subgroups, i.e. between controls and two MDMD subgroups.In contrast, nominal variable comparisons were performed using contingency table analysis, specifically the Chi-square test.In addition, Pearson's and point-biserial correlation coefficients were used to analyze the relationships between scale variables and binary variables.False discovery rate (FDR) p value correction was applied to adjust for multiplicity.The researchers used manual multiple regression analysis to investigate the impact of ACEs, NLEs, and additional demographic factors on the immune profiles.In a similar manner, these analyses examined the influence of different predictor variables, namely immune profiles, ACEs, and NLEs, on the manifestations of depression.In addition, automatic forward stepwise regressions were employed.A significance level of p = 0.05 was used to determine the inclusion of variables, and a significance level of p = 0.06 was used to determine their exclusion.This method made it easier to determine which variables should be included in the final regression model and which ones should be excluded.For each variable included in the final regression models, the standardized coefficients, t-statistics, and exact p values were calculated.In addition, we determined the F statistics, their respective p values, and the effect magnitude using the partial Eta squared.Using appropriate statistical measures, the presence of multicollinearity and collinearity in the data was extensively examined.A tolerance limit of 0.25 and a threshold for the variance inflation factor of four were used.In addition, for the purpose of this analysis, we utilized the condition index and variance proportions derived from the collinearity diagnostics table.Using the White test and a modified variation of the Breusch-Pagan test, heteroskedasticity was identified.All the www.nature.com/scientificreports/preceding analyses employed two-tailed tests.A significance level of 0.05 or lower was regarded as statistically significant.We show significance levels as exact p values, unless the p value is lower than 0.001.The latter cases are shown as p < 0.001 (or 0.000 in the partial least squares figures).As required, we utilized transformations such as logarithmic, square-root, rank-based inversed normal (RINT), and a Winsorization technique to obtain a normal distribution for our data indicators.Partial least squares (PLS) analysis was used to investigate the causal relationships between ACEs and NLEs, immunological profiles, and the physiosomatic symptoms of depression.The output variable was a latent vector derived from the diverse symptom domains, and ACE + NLEs and immune profiles as explanatory variables.In addition, the immune profiles were allowed to mediate the effects of ACE + NLEs on the clinical assessment scores.Only when both the external and internal models met the predetermined quality criteria was a comprehensive partial least squares (PLS) analysis conducted.These quality criteria are: (a) Confirmatory tetrad analysis (CTA) verifies that the latent vectors derived from the indicators have not been incorrectly specified as reflective models.(b) The blindfolding procedure reveals that the cross-validated redundancy of the construct is adequate.(c) The latent vectors exhibit strong construct and convergence validity, as indicated by composite reliability values greater than 0.80, Cronbach's alpha values greater than 0.70, and average variance extracted (AVE) values greater than 0.5.(d) At a significance level of p < 0.001, all loadings on the extracted latent vectors exceed 0.65.(e) The model fit is regarded satisfactory if the standardized root squared residual (SRMR) is less than 0.08.Consequently, a thorough pathway analysis is conducted using PLS-Structural Equation Modeling (PLS-SEM).In the analysis, SmartPLS software and 5,000 bootstrap samples were utilized.Path coefficients and their respective p values were calculated.In addition, specific indirect effects, total indirect effects (mediated effects), and total effects were calculated if the model quality data met the specified conditions.The estimated minimum sample size is 103 based on a power analysis (conducted with G*Power 3.1.9.4) using a linear multiple regression analysis with an estimated effect size of 0.111 (corresponding to approximately 10% of the variance explained), a significance level (alpha) of 0.05, and statistical power of 0.8, with three covariates.

Ethical approval and consent to participate
The College of Medical Technology at The Islamic University of Najaf, Iraq (18/2021) approved the research project.Our IRB follows the International Guideline for Human Research Safety, as well as the World Medical Association Declaration of Helsinki, The Belmont Report, the CIOMS Guideline, and the International Conference on Harmonization of Good Clinical Practice, and our study was conducted in accordance with all applicable Iraqi and international ethics and privacy laws.(ICH-GCP).

Informed consent statement
All participants and their parents or legal guardians signed a written consent form.

Socio-demographic and clinical data
Table 1 shows the sociodemographic and clinical data of the participants (controls and MDMD patients) divided into three groups using a PC extracted from 6 symptom domains (pure_depression, pure_anxiety, physiosomatic, insomnia, melancholia, and gastro-intestinal symptoms), and the major immune profiles (M1, Th-1, IRS, CIRS) and cytokines/chemokines (IL-16, TRAIL, CLC27, SCGF, M-CSF) (labeled: PC_immune + phenome) that differentiate controls from patients.We were able to extract one validated PC from these 15 variables (KMO = 0.839, Bartlett's test of sphericity = 1155.581,df = 36, p < 0.001, EV = 61.5, all loading of the 15 variables are > 0.694).Consequently, using a visual binning method the PC score was divided in three groups using − 0.606 and + 0.608.As such, the MDMD group was divided into two subgroups namely severe versus less severe MDMD and these were compared with healthy controls.Table 1 shows the mean values of the PC_immune + phenome scores.There was a strong association between these groups and the clinical diagnosis (controls versus MDMD patients): the PC-derived class with a PC score < − 0.606 were the 32 controls, while the patients were divided into those with severe symptoms and immune aberrations (labeled "severe MDMD") versus milder symptoms and immune disorders ("milder MDMD").
Table 1 shows that there were no significant differences in age, sex, education, marital status, BMI, and previous COVID-19 infection between the three study groups.Pure_depression, pure_anxiety, pure_FF, and insomnia scores were higher in patients than controls, while physiosomatic, melancholia, and GIS scores were different between the three groups and increased from controls to milder MDMD to severe MDMD.Autonomic symptoms were significantly increased in the latter group as compared to controls.There were no significant differences in the use of fluoxetine (20-40 mg per day), amitriptyline (150-200 mg per day), escitalopram (10-20 mg per day), olanzapine (5-15 mg per day), and mirtazapine (15-30 mg per day just before sleep) between the two depression study groups.There were twenty-six drug-naïve patients (never taken any antidepressant) and the medicated/ unmedicated ratio was not significantly different between both depression subgroups.Moreover, there were no significant differences in the associations between physiosomatic symptoms and immune profiles between drug-naïve patients and those using antidepressants.As a consequence, we have combined these two subgroups.

Intercorrelation matrix
Table 3 shows the correlation between the major symptom domains and the immune profiles.There were no significant associations between any of the immune profiles and autonomic symptoms.The zM1-zM2 and zIRS-zCIRS scores were not correlated with all symptom domains.The M1 classical, Th-1, IRS, CIRS, and Immu-NT profiles were significantly correlated with all symptom domains (except autonomic symptoms).The alternative M2 profile was associated with all symptom domains, except pure FF, GIS, and autonomic symptoms.Th-2 was significantly correlated with pure_depression, pure_anxiety, melancholia, and insomnia.Th-17 was associated with all symptom domains except pure_FF and autonomic symptoms.The zTh-1-zTh-2 profile was significantly correlated with all domains, except melancholia, insomnia, and autonomic symptoms.

Results of multiple regression analysis
Table 4 shows the results of different multiple regression analyses with the symptom domains as dependent variables and the cytokines/chemokines/growth factors as dependent variables, while allowing for the effects of age, sex, smoking, education, and BMI.Since we found that multiple CIRS cytokines were significant in the final regression analysis we also entered the CIRS score as an additional explanatory variable.Regression #1 shows that 48.9% of the variance in the pure_depression score was explained by IL-16 and SCGF (both positively) and CCL3 (inversely).Regression #2 shows that 45.4% of the variance in the pure_anxiety score was explained by IL-16 and IL-17 (both positively) and CCL3 (inversely).Regression #3 shows that 59.0% of the variance in physiosomatic symptoms was explained by IL-16, IL-8, and male sex (all positively associated) and CCL3 and sIL-1RA (both inversely).Figures 1 and 2 show the partial regression of the physiosomatic symptom domain on IL-16 and CCL3, respectively.We found that 46.5% of the variance in the pure_FF score was explained by IL-16 and TRAIL (positively) and inversely by the CIRS score (regression #4).We found that 53.2% of the variance in www.nature.com/scientificreports/melancholia was explained by the regression on IL-16 and SCGF (both positively) and CCL3 (inversely).A large part of the variance in GIS was predicted by IL-16, CCL1, TRAIL, and IL-8 (all positively) and CCL3 (inversely).Up to 54.1% of the variance in the insomnia score was predicted by IL-16 and SCGF (both positively) and sIL-1RA (inversely).Only a small part of the variance in autonomic symptoms was explained by TRAIL.

Effects of ACEs and NLEs on the phenome are mediated by immune aberrations.
Since there are associations between ACE + NLEs and the immune profiles as well as symptom domains 33 , we also examined whether the regression of the symptoms domains on the immune products was influenced by ACE + NLEs.Table 5 shows the results of multiple regression analyses which included the significant effects of ACE + NLEs on the phenome data.Thus, the ACE + NLEs score had a significant effect on all symptom domains, except autonomic symptoms.The results of this Table also show that cytokines/chemokines/growth factors had significant effects on the symptom domains beyond the effects of ACE + NLEs.In all cases, the explained variances were significantly enhanced by entering the ACE + NLEs values with around 3% (variance explained), except for GIS (no significant improvement in explained variance).

Results of PLS analysis
The first PLS model conceptualized the phenome of depression as a latent vector extracted from physiosomatic and FF symptoms, as well as pure_depression, pure_anxiety, GIS symptoms, melancholia, and insomnia.The ACEs + NLEs index was composed as a formative model.Based on the results of the multiple regression and correlation analyses, we entered the most significant cytokines/chemokines/growth factors predicting the phenome, and in addition included the effects of ACE + NLEs on the immune variables (together with age and sex).The  www.nature.com/scientificreports/composite reliability (0.956), Cronbach's alpha (0.945) and AVE (0.759) of the phenome latent vector were more than adequate, as was the model quality fit with an SRMR of 0.047.All Q 2 predicted values for the manifest and latent variables were positive using PLSpredict, indicating that the model outperforms the naivest benchmark as shown in Fig. 3.We observed that the phenome of depression was best predicted by ACE + NLEs, IL-16, TRAIL (all positively), and sIL-1RA and CCL3 (both inversely).The latter five predictors explained 69.6% of the variance in the phenome of depression.In addition, part of the variance in IL-16 (26.8%),TRAIL (19.4%) and sIL-1RA (6.4%) was explained by ACE + NLEs.Age was associated with CCL3 and sex with TRAIL.As such, there were significant specific indirect effects of ACE + NLEs on the phenome which were partly mediated by IL-16 (t = 4.40, p < 0.001) and TRAIL (t = 2.16, p = 0.040).The direct and total indirect (t = 5.99, p < 0.001) effects of ACE + NLEs on the phenome explain the strong total effects of ACE + NLEs on the phenome (t = 13.21,p < 0.001).Consequently, we have constructed a second PLS model by integrating different cytokines/chemokine/growth factors into one latent vector, labeled "immune network, " including IL-16, TRAIL, CCL27, M-CSF, and SCGF.This factor showed adequate convergence and construct validity as indicated by an AVE value of 0.736, composite reliability of 0.933, and Cronbach's alpha of 0.910.Moreover, the overall model fit was adequate with an SRMR value of 0.043.We found that 68.8% of the variance in the phenome of depression was explained by this immune network (positively associated) and CCL3 (inversely associated).ACE + NLEs had significant effects on the phenome, comprising direct effects (see Fig. 4) as well as specific indirect effects (t = 6.55, p < 0.001) which Table 4. Results of multiple regression analysis with symptom domains as dependent variables and cytokines/ chemokines/growth factors as explanatory variables.IL interleukin, SCGF stem cell growth factor, sIL-1RA soluble IL-1 receptor antagonist, sex males = 1, females = 0, CIRS compensatory immunoregulatory system, TRAIL tumor necrosis factor-related apoptosis-inducing ligand.www.nature.com/scientificreports/were partially mediated by the immune network.Both effects of ACE + NLEs culminated in strong total effects (t = 13.41,p < 0.001).On the other hand, there was no significant effect of age on the phenome (t = 1.62, p = 0.106).

Discussion
The

physio-affective phenome of depression
The first major discovery of this study indicates that within MDMD patients, the physiosomatic, FF, and GIS symptom domains are associated with the same factor as pure depression, pure anxiety, melancholia symptoms, and insomnia.This suggests that the various symptom domains are closely interconnected expressions of the same latent construct, specifically the phenome of depression.Previous studies have demonstrated that physiosomatic symptoms are associated with the same phenome factor as pure depressive and anxiety symptoms 7,46 .
To emphasize the influence of interconnected affective and physiosomatic symptoms on the manifestation of MDD, this construct has been renamed as the physio-affective phenome 7,47 .Furthermore, various medical and Table 5. Results of multiple regression analysis with symptom domains as dependent variables and cytokines/ chemokines/growth factors, adverse childhood experiences (ACEs) and negative life events (NLEs) as explanatory variables.IL interleukin, SCGF stem cell growth factor, sIL-1RA soluble IL-1 receptor antagonist, sex males = 1, females = 0, CIRS compensatory immunoregulatory system, TRAIL tumor necrosis factorrelated apoptosis-inducing ligand.neuro-psychiatric disorders exhibit a latent construct encompassing symptoms of depression, anxiety, and physiosomatic manifestations.This phenomenon can be observed in conditions such as acute COVID-19 infection 48 , Long COVID 49 , relapsing-remitting multiple sclerosis (RRMS) 50 , schizophrenia 51 , and epilepsy 52 .Furthermore, FF symptoms belong to the same latent construct which loads highly on affective and physiosomatic symptoms.The aforementioned interconnected relationships among FF symptoms and the phenome of depression were previously identified in the context of MDD 7 , acute COVID-19 infection 53 , RRMS 50 , Long COVID 54 , schizophrenia 55 , and end-stage renal disease 56 .
The interrelationship between affective, physiosomatic, and FF symptoms is well-established, indicating that these symptoms are interconnected and should not be considered distinct phenomena.This implies that major depressive disorder (MDD), bipolar disorder, and various psychiatric disorders exhibit a physio-affective phenome.These findings also suggest that the physiosomatic, FF and GIS symptom domains, as well as pure depression and anxiety, melancholia, and insomnia, are associated with similar underlying mechanisms that involve both central and peripheral aberrations.This observation suggests that MDMD is a condition that affects the entire body and that affective symptoms are merely one component of the comprehensive clinical presentation.

The physio-affective phenome is predicted by immune activation
The results of the present study show that physiosomatic and FF symptoms of MDMD are best predicted by a) IRS, Th-1 activation, increased Immu-NT, Th-1 polarization, and M1 activation; b) a combination of increased IRS products (especially IL-16, IL-8 or TRAIL), lowered CIRS components (sIL-1RA), and lowered levels of a chemokine (CCL3), and c) an immune network comprising IL-16, TRAIL, CCL27, M-CSF and SCGF (all positively) and CCL3 (inversely).These results indicate that T cell activation (IL-16, Th-1 polarization) is a significant factor in the development of physiosomatic, FF, and GIS symptoms.
IL-16 is elevated in major depressive disorder 57 and is a crucial cytokine in FE-MDMD 24 .IL-16 and other proinflammatory cytokines are associated with elevated depression and neuroticism scores in hepatitis C patients receiving cytokine-based immunotherapy 58 . 59,60.IL-16 induces T cell signaling via the CD4 molecule on Th cells, resulting in the upregulation of activation markers such as HLA-DR + and CD25 + cells.Prenatal depression is associated with elevated IL-16 levels in the prefrontal cortex and hippocampus in animal models 61 and contributes to neuroinflammation in animal models of autoimmune encephalomyelitis by activating T cells 62 .CCL27 or CTACK (cutaneous T cell-attracting chemokine) is a pro-inflammatory chemokine that attracts memory T cells (T lymphocyte-associated antigen) and this chemokine is essential for T cell-mediated cutaneous inflammation 63,64 .Such interactions may explain the co-occurrence of depression and "psychocutaneous disorders" like psoriasis, which is accompanied by elevated CCL27 levels 65 .Moreover, CCL27 is a key regulator of immune homeostasis in mucosal tissues and the epidermis 66 , and CCL27 KO-mice exhibit CR10 + T cell infiltration in the reproductive tract and lungs 66 .Another cytokine that is strongly associated with physiosomatic and FF symptoms is TRAIL.TRAIL signals apoptosis by binding to death receptors and activating the caspase pathway 67 .This cytokine has strong neurotoxic effects, even leading to neurodegeneration and neuronal death, mediated by the TRAIL-R2/DR5 receptor and caspase activation 68 .Recent evidence suggests, however, that TRAIL may modulate the IRS and T cell activation and shows anti-inflammatory effects 67 .TRAIL inhibits the activation of T cells, indicating that this cytokine functions as an activator of an anti-inflammatory pathway independent of apoptosis.
Many patients with MDD exhibit T cell activation, as measured by flow cytometry and T cell activation markers including HLA-DR+ and CD7+ CD25+ T cells 69,70 .We have argued, based on these and other findings, that cell-mediated activation, and specifically T cell activation, is a key factor in MDD 37 .CD3+ CD71+, CD3+ HLADR+, and CD4+ CD71+ T cell phenotypes have been demonstrated to differentiate MDMD from SDMD 71 .In addition, the numbers of CD3+ CD40L, CD4+ CD40L, CD4+ HLADR+, and CD8+ CD40L+ cells are substantially higher in MDMD than in healthy controls, whereas SDMD patients occupy an intermediate position 71 .
Importantly, physiosomatic and FF symptoms are strongly associated with Th-1 cytokine profiles, but not with Th-2 profiles.Th-2 cytokines (such as IL-4) have negative immune-regulatory and anti-inflammatory effects and are essential for CIRS activity.As a result, a relative decrease in Th-2 functions may result in attenuated immunoregulatory Th-2 effects, thereby enhancing the deleterious effects of Th-1 and IRS activation.Similarly, our results demonstrate that decreased CIRS activity plays a crucial role in FF symptoms.Even more significant is the fact that relative decreases in sIL-1RA levels (another component of CIRS) contribute to physiosomatic symptoms and the physio-affective phenome of depression.There is evidence that both MDD and MDMD are associated with elevated sIL-1RA serum/plasma concentrations 24,72 .Maes and Carvalho (2018) report that soluble sIL-RA may inhibit both IL-1α and IL-1β signaling, thereby exerting potent anti-inflammatory effects.Recent findings indicate that the phenome of MDMD is significantly predicted by activated T and B cells, whereas T regulatory cells (such as CD25 + FoxP3 + GARP +) are inversely associated with the phenome 71 .Based on the IRS/CIRS theory of depression 73 , we may conclude that increased Th-1 (IL-16), TRAIL and M1 activities, and relative deficiencies in CIRS (including Th-2 cytokines and sIL-1RA) are major contributors to physiosomatic and FF symptoms.
The analyses that we conducted in the current study suggest that a decrease in CCL3 may contribute to the severity of physiosomatic and FF symptoms, which is another potentially significant finding.CCL3 (or macrophage inflammatory protein) is a pro-inflammatory chemokine that recruits activated immune cells to inflammatory sites and increases the expression of pro-inflammatory M1 cytokines including IL-1β, TNF-α, www.nature.com/scientificreports/and IL-6 74 .Previously, we found that whole blood-stimulated CCL3 production in the culture supernatant of patients with MDMD was significantly lower than that of healthy controls 75 .In addition, CCL3 was one of the few cytokines/chemokines whose serum concentrations were lower in the serum of MDMD patients compared with controls 24 .Consequently, we can hypothesize that relative decreases in CCL3 secretion may play a role in the depression phenome.CCL3 mediates fever independent of prostaglandin (PGE2) and elevates PGE2 in cerebrospinal fluid [76][77][78] .Endogenous PGE2 possesses context-dependent anti-inflammatory properties and contributes to the resolution of inflammation 79,80 .In addition, CCL3 has neuroprotective properties, such as protecting dopaminergic neurons, and influencing traumatic brain injuries 81 .

Effects of ACE and NLEs on the phenome are mediated by the immune network
The third major finding is that the combined independent effects of immune aberrations and ACE + NLEs significantly predict the physiosomatic and FF symptoms of MDMD, and the effects of ACE + NLEs on the phenome are partially mediated by increased IL-16 and TRAIL, whereas ACE-activation of sIL-1RA levels may have protective effects.Our second PLS analysis revealed that ACE + NLEs exert a substantial influence on an immune network composed of IL-16, CCL27, TRAIL, MCSF, and SCGF, which partially mediates the effects of ACE + NLEs on the phenome, including physiosomatic and FF symptoms.SCGF (CLECC11a) activates different hematopoietic progenitor cells and activates cell signaling networks including glycogen synthase kinase 3, β-catenin and Wnt pathways, which play a major role in MDD, physiosomatic and FF symptoms 2,82 .M-CSF is a macrophage differentiation factor that is essential for microglial activation and neuroinflammation 83,84 .Together, ACEs and NLEs appear to induce a pro-inflammatory network composed of cytokines, chemokines, and growth factors.These findings extend those of a previous paper 31 , which demonstrated that ACEs elicit a network of cytokines/chemokines/growth factors measured in the culture supernatant of stimulated whole blood from MDMD patients.In fact, both the present investigation (conducted on Iraqi patients) and our previous report (conducted on Thai patients) 2 demonstrate that ACEs activate M1, Th-1, Th-17, IRS, and CIRS profiles and significantly increase Immu-NT.In addition, ACEs were associated with elevated levels of sIL-1RA, IL-9, IL-12p70, PDGF, and TNF-α in both studies.
Maes et al. 2 demonstrated that ACEs induce stimulated but not unstimulated levels of cytokines/chemokines/ growth factors, suggesting that the latter factors are sensitized in MDMD patients.Future stressors could then reactivate these sensitized networks, according to the hypothesis.Since NLEs increase the serum levels of IL-16, CXCL12, M-CSF, SCGF, TRAIL, IL-18, PDGF, sIL-2R, sIL-1RA, SCF, and IL-9, and ACE + NLEs have a greater effect on CXCL12 and M-CSF than ACEs, we can conclude that NLEs may indeed reactivate the sensitized ACE-induced immune responses.However, ACE + NLEs also had significant direct effects on the MDMD phenome, indicating the involvement of other pathways, such as gut-brain pathways, gut-dysbiosis, oxidative stress pathways, and decreased antioxidants and neurotrophic protection 4,16,29,71 .

Limitations
This paper might have been more intriguing had we assessed additional significant biomarkers of MDD, such as oxidative and nitrosative stress, which are linked to IRS activation 1 .An additional concern pertains to the utilization of a machine-learning-derived classification for MDMD, a more severe subtype of the disorder.Consequently, the generalizability of the results from this investigation to moderate subtype of depression (SDMD) or the overarching group MDD 1,12,18 may be limited.Moreover, it is imperative that our findings be replicated in diverse nations and cultures.

Conclusions
The findings of this study indicate that physiosomatic, FF and GIS symptoms are strongly associated with depression, anxiety, melancholia, and insomnia domains of MDMD.The first factor derived from these six different domains is designated as the physio-affective phenome of depression.A larger part of the variances in the physioaffective phenome, and the physiosomatic FF and GIS symptom domains is explained by immune variables.Specifically, IL-16, TRAIL and IL-8 have a positive correlation with these symptoms, while CCL3 and sIL-1RA exhibit an inverse correlation.The application of PLS analysis reveals that the combination of ACEs and NLEs has a significant impact on the physio-affective phenome.This influence is partially mediated by an immune network consisting of interleukin-16, CCL27, TRAIL, M-CSF, and SCGF.Overall, the physiosomatic and FF symptoms of FE-MDMD are attributed, in part, to Immu-NT linked to Th-1 polarization and M1 macrophage activation, as well as a decrease in CIRS protection.IL-16, Th-1 polarization, M1 activation and Immu-NT are new drug targets to treat the physiosomatic and FF symptoms of MDMD, and to desensitize the upregulated cytokine network and deprogram the detrimental effects of ACEs and NLEs. https://doi.org/10.1038/s41598-024-57350-1

Table 3 .
Intercorrelation matrix between immune profiles and symptoms domain scores.